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The derivation of a Moving Boundary Approximation or of the re- 
sponse of a coherent structure like a front, vortex or pulse to external 
forces and noise, is generally valid under two conditions: the existence 
of a separation of time scales of the dynamics on the inner and outer 
scale and the existence and convergence of solvability type integrals. 
We point out that these conditions are not satisfied for pulled fronts 
propagating into an unstable state: their relaxation on the inner scale 
is power law like and in conjunction with this, solvability integrals di- 
verge. The physical origin of this is traced to the fact that the important 
dynamics of pulled fronts occurs in the leading edge of the front rather 
than in the nonlinear internal front region itself. As recent work on the 
relaxation and stochastic behavior of pulled fronts suggests, when such 
fronts are coupled to other fields or to noise, the dynamical behavior 
is often qualitatively different from the standard case in which fronts 
between two (meta)stable states or pushed fronts propagating into an 
unstable state are considered. 
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1 Introduction 



For a pattern in two or more dimensions that naturally can be divided into do- 
mains and "domain walls" separating them, a much used analytical approach 
is a moving boundary or effective interface approximation |], ||, ||, |], |(| [?J |§. 
This seems appropriate, when the width of the domain wall, front, interface, 
or transition zone is much smaller than the typical length scale of the pattern 
and when the dynamics of the pattern on long space and time scales occurs 
through the motion of these interfaces. The moving boundary approximation 
amounts to treating these fronts or transition zones as a mathematically sharp 
interface or boundary. In other words, their width is taken to be zero and 
their internal degrees of freedom are eliminated. We shall henceforth use the 
word boundary or interface to denote this zero width limit and use the word 
front when we look at a scale where its internal structure can be resolved. 

Moving boundary approximations (MBA's) are ubiquitous in the theory 
of pattern formation: they arise in most analytical approaches to late stage 
coarsening ||, 1Q[ , in the analysis of interface dynamics in dendritic growth 



and viscous fingering [II], [T^, [13], [TJ|, [T^, [IB], [T7|] , step dynamics at surfaces [fTS], 



19| , |20l, thermal plumes [[H], 22]] , in chemical wave dynamics [p3fl , combustion 



fronts 0, etc. 

The main physical idea underlying the derivation of a MBA is that the 
front itself can on large length and time scales be viewed as a well-defined 
coherent structure which can be characterized by its coordinates and a few 
effective parameters, such as its velocity or a mobility coefficient. This idea 
plays a role for many coherent structures, like vortices, or pulse-type solutions 



like sources, sinks, solitons, etcetera |24 , 25] , |2q , 27 , |2q] . The response of a co- 



herent structure to an external driving force or noise [^9], |30|] or the interaction 
between them can frequently be derived by a perturbative expansion about 
the isolated coherent structure solution. Often the effective parameters (a dif- 
fusion coefficient, a mobility or an effective interaction force) can be derived 
from a solvability condition. A solvability condition expresses that a linear 
equation of the form L <pi — gi, where the linear operator L results from lin- 
earizing about the isolated coherent structure solution, is solvable provided g\ 
is orthogonal to the kernel (null space) of L. In other words, the requirement 
for such an equation to be solvable is that that gi is orthogonal to the left zero 
mode x of L. 

Although this is hardly ever mentioned explicitly, there are two important 
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implicit assumptions underlying such approximations, namely (a) that there 
is a separation of time scales between the motion of the front as a whole 
and its internal dynamics, and (b) that the internal dynamics of the front is 
determined by the nonlinear front region itself, so that the solvability type 
integrals are dominated by the contributions from this finite region, and hence 
donot diverge. 

The issue that we address in this paper is that while the above conditions 
are satisfied for the familiar MBA for bistable fronts and also for so-called 
pushed fronts, they are not for so-called pulled fronts. We will indeed discuss 
several related properties of pulled fronts which bear on this: (i) the divergence 
of the solvability integrals, with the concommittant breakdown of a MBA or of 
the derivation of the response functions of the front, like a diffusion or mobility 
coefficient; (ii) the shift of the dynamically dominant zone from the interior to 
the leading edge of the front, that causes the solvability integrals to diverge; 
(Hi) the fact that the stability spectrum of planar pulled front solutions is 
gapless; (iv) the recently discovered universal slow power law relaxation of 
planar pulled fronts . We will initially focus our discussion on the derivation 
of a MBA, but as we shall see our observations and conclusions apply equally 
well to essentially any perturbative analysis of a pulled front. 

The crucial feature of the standard moving boundary problem is that the 
boundary conditions are local in space and time — e.g., the growth velocity of 
an interface is a function of the instantaneous local temperature and curvature 
of the interface. Usually, some of the boundary conditions are associated with 
conservation laws, like the conservation of heat, and so they can often be 
guessed from physical considerations. 

If there is a separation of spatial scales, then such a MBA applies only if 
there is also a separation of time scales between the internal dynamics of the 
front and the dynamics of the outer bulk fields. E.g., if the internal front modes 
relax on a time scale r, and one considers a front of width W, propagation 
velocity v and typical curvature k, then a MBA becomes appropriate in the 
regime kW <C 1, vkt -C 1. Such a well-defined relaxation time r of a front on 
the inner scale actually exists only if the relaxation is exponential in time. In 
this case, r is the inverse of the gap in the spectrum of the stability modes of the 
planar front. Just like multiple scale and amplitude expansions [32, 33, [34j are 



based on projecting all rapidly decaying gapped modes onto the slow one (the 
center manifold), the MBA or effective interface approximation can be thought 
of as projecting a problem with fronts onto the slow interfacial dynamics. 
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However, if the stability spectrum of the planar front is gapless, the internal 
modes of the front relax algebraically in time. Thus there is no characteristic 
time t for the internal modes, no separation of time scales and no standard 
MBA, no matter how thin the front is. The internal dynamics of such a pulled 
front is actually slaved to the evolution of its leading edge on the outer scale, 
which motivates the term "pulling". Note that despite its different temporal 
behavior, it is not at all visible on an instantaneous picture of a front, whether 
it is bistable, pushed or pulled. 

In a problem where the starting equations are partial differential equations, 
the derivation of the MBA can often be done analytically using by now stan- 
dard methods. One should keep in mind, however, that MBA's can be equally 
powerful in situations where the approximation can not be derived cleanly by 
starting from a partial differential equation and applying standard methods. 
E.g., in crystal growth the interfacial boundary conditions are determined on 
a molecular scale, where for a rough interface the molecular processes are so 
fast that after some coarse graining, we can describe the interface for many 
purposes as a sharp interface whose response to changes in temperature and 



concentration are instantaneous |3q| . Similar considerations apply to coarsen- 
ing interfaces or combustion fronts. 

In the next section, we will first summarize the necessary essentials of the 
stability and relaxation properties of pulled fronts. Then, in section 3 we 
illustrate the issue by following the standard derivation of a MBA for the 
type of coupled equations that have in recent years been used in a phase- 
field type formulation of solidification problems. In section 4 we then discuss 
the conditions under which such a type of analysis applies in more detail, to 
identify the difficulties that arise when the front dynamics on the inner scale is 
changed from the usual bistable or pushed case to pulled. We then in section 
5 generalize our findings to equations with higher derivatives and to coupled 
equations, that create uniformly translating fronts. We show that the usual 
route of deriving solvability conditions does work in general for bistable and 
pushed fronts, but not for pulled fronts. 
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2 Pulled fronts: Properties and statement of 
the problem 

When one considers a linearly unstable state, even a small perturbation about 
this unstable state grows out and spreads. We will confine our analysis to fronts 
emerging from a localized initial perturbation of the unstable state. One can 
calculate the asymptotic linear spreading velocity v* of such a perturbation 
simply from the linear dispersion relation u(k) of the unstable modes according 

to [06i m 



du(k) 



dk 



lmuj(k* 



Im k* 



We furthermore will confine ourselves in this paper to fronts which asymptot- 
ically are uniformly translating. For these, uj* and k* are purely imaginary, 
and we use the notation k* = iX*. If the above equation admits more than 
one solution, the one corresponding to the largest value of v* is the relevant 
one. We refer for the derivation of these results to our recent paper |3l|, which 



we will quote as paper I below. Pulled fronts are those for which the asymp- 
totic spreading velocity v as of the nonlinear front equals this linear spreading 
velocity v*\ v as = v* f 



39], m [41], [31| . A number of model equations for 



which fronts are pulled are discussed in paper I, but they also arise in the 
analysis of more complicated situations like pearling [Q, the Couette- Taylor 
instability |43[| , Rayleigh-Benard convection JPJ, the instabilities of wakes of 
bluff bodies, leading e.g. to von Karman instabilities the emergence of 
global modes H^j, liquid crystals fT7| , streamer discharge patterns the 
competion of domains in the Kupers-Lortz instabilility fi9 |, the emergence 
of domains near structural phase transitions [pTj], polymer patterns su- 



propagation 



perconducting fronts p^fl , error propagation |53[, deposition models p3fl , step 
chaotic fronts in the complex Ginzburg-Landau equation 

and the 



18 



57 1, renormalization group analysis of disorder models |58 



analysis of the Lyapunov exponents in kinetic models |59| . 

Fronts which propagate into an unstable state always are pulled if all the 
nonlinearities suppress the growth. If not all of them do, the asymptotic fronts 
speed v as may become larger than v*: v as = > v*. The relaxation of such 



"pushed^ fronts |38], fy]> |3_TJ] is exponential with a characteristic relaxation time 
r, that is finite KD, [31]. As we will discuss, for these the same perturbative 



schemes apply as for the familiar bistable fronts, and likewise for these a 
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standard type MBA can be derived. 

In paper I, we have shown that when a pulled front grows out of sufficiently 
steep initial conditions (decaying into the unstable state at least as e~ Xx for 
x — > oo with some A > A*), then the velocity of a front obeys a universal 
power law relaxation given by 

v(t) = V* + X(t) , (2) 



where 



D _<dM*) 



2dP 



(4) 



is real and positive for uniformly translating front solutions. For the front 
profile, a similar power law relaxation holds, and the extension of these results 
to one-dimensional pattern forming fronts is given in p7[] . The analysis reveals, 
that the power law relaxation emerges from the dynamics of the foremost part 
of the front where the dynamics is governed by the equations linearized about 
the unstable state. The dynamics in the nonlinear region is essentially slaved 
to this so-called leading edge. The very slow 1/t power law relaxation of pulled 
fronts without characteristic time scale obviously implies that the separation 
of time scales, which is necessary for a MBA to be applicable, is missing. 

While from this perspective it is already intuively obvious that a standard 
perturbation theory or MBA does not apply to pulled fronts, the arguments 
underlying the separation of time scales are hardly ever discussed explicitly in 
the literature on the derivation of a MBA. The purpose of this article there- 
fore is to point out where the standard derivation breaks down and how this 
emerges at a more formal level. In such an approach, one generally encounters 
solvability type integrals of the form 



or generalizations thereof — see e.g. (p2|) , (|27|) or (p7|). Here £ = x — vt is a 
frame moving with the front with speed v, and $o(0 is the associated planar 
front solution. The translation mode (9^$o is a right zero mode of the linear 
operator L emerging from linearization about the asymptotic front Qq, and 
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e v £ c^$ is a left zero mode of this operator. As we shall see, such solvability 
type integrals are well-defined and finite for bistable and pushed fronts, but 
diverge for pulled fronts, since the integrand does not converge for £ — > oo. 
In a way, the solvability integral still correctly distributes its weight over the 
dynamically important region, but for a pulled front, this region becomes 
semi-infinite, and therefore the integral diverges. Our discussion also shows 
why introduction of an ad-hoc cutoff in these integrals — an approach that 
has sometimes been considered in the literature — does not necessarily cure 
the problem. 



3 The derivation of a MBA from a phase field 
model 

In this section, we first follow the standard derivation of a moving boundary 
approximation (MBA) from a phase field model to highlight the assumptions 
and approximations along the way. We then analyze why and how the approx- 
imation breaks down for pulled fronts. 

As an example, we study the "phase field model" 

Ou/dt = aV 2 u + d(j)/dt , (6) 
ed<j)/dt = e 2 V 2 (f) + f(<f>,u) , (7) 
where f(</>, u) = <f) (1 - <p) - \u + <p) , A > . (8) 

In the limit of zero front width e — > 0, this model for appropriate parameters 
a, fi and A reduces to a moving boundary approximation for a solidification 
front, where we can think of <fi as the order parameter field, while u plays 
the role of the temperature. then varies from the stationary "liquid-like" 
solution <p in one domain to another "solid-like" solution « 1 in the 
other domain. Note that in contrast to d<p/dt in (|7|) has a coefficient e, 
not e 2 . This allows the front to have a velocity of order unity, so the velocity 
is nonvanishing already in the lowest order perturbation theory 0(e°). The 
d(j)/dt on the r.h.s. in (^) models the generation of latent heat in the interfacial 
zone where <fi changes rapidly. 

Other choices for /(</>, u) can be found in the literature [0, £|, || but the 



form (pf) is most convenient for our present purpose. / can be considered as 
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the derivative of a "free energy" F, 

/OM = q^- > ( 9 ) 

- Xu) 4> 2 Xu) 3 4 

r \(b,u) = . 

W j 2 3 4 

Since w varies on spatial and temporal scales of order unity, let us treat it as 
a constant for a moment on the small length scale e on which varies, and 
let us define p = fi — Xu. The connection with the phase field models for 
solidification is closest in the range — 1 < p < 0, when the function F (which 
is like a Ginzburg-Landau free energy density), has two minima at = and 
at = 1. When p = -1/2, then F(0,u) = F(l,u) = and the two "phases" 
= and = 1 are in equilibrium. So if we choose the bare parameter 
\x = —1/2, then u = corresponds to the melting temperature, where (|?|) 
admits stationary front solutions with velocity v — 0. For fi = — 1/2 but u 
nonzero, the minima of / shift relative to each other, and the order parameter 
front ([7D moves. When u is positive, the liquid like minimum at = is the 
absolute minimum of F, and for u negative the solid like minimum at = 1 is 
the absolute one. The front then will move such that the state with the lowest 
free energy extends. For p, > the state = is linearly unstable; so we 
then deal with fronts propagating into unstable states which are pushed for 
< p < 1/2 and pulled for p > 1/2 |3T|. Though the interpretation 



of the model as a solidification model might be lost, we will illustrate the 
derivation of a MBA as a function of p for this example, and we will find 
that the method breaks down at the transition from pushed to pulled fronts 
at p = 1/2. 

Let us now trace the steps of the approximation in more detail. The field u 
(||) varies on a spatial scale of order unity, and the field (0) on a spatial scale 
of order £< 1. A moving boundary approximation consists of first matching 
an inner expansion of the problem on scale e to an outer problem on scale 
1, and then letting e — > such that an effective moving boundary problem 
on the outer scale results. In the limit of e — >• 0, the interface might have a 
nonvanishing velocity and curvature on the outer length scale, so we allow for 
v = 0(e°) andK = 0(e°) @. 

Let us for simplicity consider the problem in two spatial dimensions (x,y). 
On the outer scale, the fields are expanded in powers of e as 

u(x,y,t) = u {x,y,t)+eu 1 (x,y,t) + ... , (10) 
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4>(x,y,t) = 4> (x,y,t)+£(f) 1 (x,y,t) + ... . 



For a further analysis of these equations on the outer scale and their matching 
to the inner scale, we refer to the literature. Here we focus on the analysis 
of the 0-front on the inner scale. First a coordinate system moving with the 
front is introduced, where s measures the arc length of the interface in the 
tangential direction, and £ the direction in which <p varies and propagates. We 
put, e.g., £ = at the place where <ft — 1/2- The coordinate £ in the direction 
normal to the front is scaled with a factor e, since the front width will be of 
order e in the limit e — > 0. However, the coordinate s is not scaled: along the 
front, the variation is assumed to be simple on length scales of the order of 
unity. For the inner expansion of the fields, one then writes 

u(x,y,t) = U {^s,t)+eU 1 ^,s,t) + ... , (12) 
4>{x,y,t) = %^,s,t)+e^,s,t) + ... . (13) 

The choice of coordinates can be illustrated when we consider a weakly curved 
front which locally propagates with a velocity v(s, t) in the x direction, so that 

x-X(s,t) 



y 



X(s,t) = x + J dt' v{s,t') 



(14) 



In general, the front is curved and has a velocity v and curvature k which varies 
locally but on the outer time scale t and spatial scale s. They are therefore 
are expanded as 

v(s,t) = v (s,t) +e vi(s,t) + . . . , (15) 
K(s,t) = Ko(s,t) +e Kt(s,t) + . . . . (16) 

The differential operators in (0) then in the interior coordinates (£, s) have the 
e expansion 



d 



e 



2 V 2 





"Ft 



i d 

v + e Vl + ... — + 0(e 2 ) 



d 2 d 



(17) 



(18) 



Inserting the expanded operators into (0) and ordering in powers of e yields 
in order e° 

— % + v — % + /($ , U o )=0, (19) 
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where Uq is essentially constant on the inner scale £. In order e 1 one finds 



n 9 * 5 * df($ ,U) 



U-i 



Uo 



with the linear operator 



/ d 2 d df(*,U ) 



(20) 



(21) 



its 



Note that since $o is a solution of the ordinary differential equation (|1 
time dependence occurs solely through the variation the U field. 

Equation ( p0|) is an inhomogeneous linear differential equation for the un- 
known field <3>i. If one has a left zero mode x(0 °f L such that xL = = L^x 
where is the adjoint operator defined through partial integrations, then ( PD[ 
can be evaluated with a so-called "solvability" analysis by projection onto %: 



<9£ 



d£ x 



df(%,U) 



dU 



U 



d£ X 



dt 



(22) 



There are clearly two important conditions for the identification of (22) with 
the common solvability condition: If the scalar products with x exist, and if 
the temporal derivative d t &o of the zero order solution (p~9|) can be neglected, 
then (|22|) expresses the first order velocity correction vx as a function of the 
local curvature k , of the outer temperature field duf Ux and of the zero order 
solution $ (0). It is exactly at these two points that the analysis breaks 
down for pulled fronts. The violation of these conditions always happens 
concomitantly, as they are physically related. 

Let us construct the left zero mode x explicitly: It is well known, that the 
right zero mode of L is the mode of infinitesimal translation d^o'- L d^o — 0. 
Since L is nonhermitian, the left zero mode of L is a right zero mode of the 
adjoint of L, 



£ f x(0 = o 



(23) 



and x <V^o(0- However, the left zero mode x can be obtained by noting 
that the transformation 



L<p = L(j> 



(24) 
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with 

/ ft 2 Rf($> TIA v 2 \ 

(25) 



4 

$0 



d£ 2 <9$ 

turns the problem into a hermitian eigenvalue problem. As a result the left 
zero eigenmode x °f L is equal to the right zero eigenmode e v ^ 2 d^o of L. 
Transforming back to L, this yields for the left zero mode 

X = e* d^ , (26) 

as can also be verified by substitution. If we may ignore the term associated 



with the time derivative <9t$o and insert the expression for x into (22) we find 

I-^e< (^) 

If we furthermore ignore the term due to the coupling to the u field, the 
expression v\ = —k,q is the familiar result of motion by mean curvature first 



derived within the context of continuum models by Allen and Cahn |6(], [10| . 

The structure of the solvability analysis is generic for the perturbative ex- 
pansion about a uniformly translating front. Although we have only considered 
the simplest type of model, and although refinements are possible 0, Eq. fl27|) 
captures the basic structure of the expression that one obtains in lowest order 
in a MBA: the relations between the velocity, curvature and temperature field 
u of the front, which play the role of boundary conditions for the outer fields 
at the boundary in the zero width limit e — > 0, contain solvability integrals 
of the form / d£ e v ^ (<9g$ ) 2 - (Note that d v f in (^?p contains a factor $ , 
that for £ — > oo decays essentially like d^o.) Solvability integrals of this type 
essentially arise in any type of perturbative calculation, since they just express 
the solvability condition of the linear perturbation problem = gi(£) : 

the inhomogeneous term <?i(£) has to be orthogonal to the left zero mode x of 
the linear operator L. 



4 Violation of the two conditions underlying 
the MBA for pulled fronts 

We now discuss the conditions under which the MBA can be derived along the 
lines sketched above in more detail. 
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Consider first the condition concerning the separation of time scales. For 
fronts between two linearly stable states (with — 1 < Jjl < in /), there is such 
a separation between the inner dynamics of the front and its displacement: In 
this case, the stability spectrum of planar front modes has a gap [35| , and all 
internal eigenmodes decay as e~ Wn */ £ with eigenvalues io n > u)\ = 1/r = 0(1). 
Thus, in the limit e — > 0, there is a clear separation of inner and outer time 
scales, and the adiabatic approximation (^7|) is justified on the outer time 
scale of order unity. Moreover, as discussed in paper I, for pushed fronts 
propagating into an unstable state the stability spectrum is also gapped, and 
therefore the separation of timescales necessary for the MBA to apply, does 
hold. However, the stability spectrum of pulled fronts is gapless, and as Eq. 
(0) of section II illustrates, pulled fronts show indeed a power law convergence 
to their asymptotic speed v*. Clearly, then, the standard derivation of a MBA 
does not apply to pulled fronts. 

The same conclusion also emerges from the properties of the solvability 
integrals themselves. For increasing v, the exponential factor e v * enhances the 
value of the integrand for large positive £, while suppressing the integrand for 
large negative £. We therefore now turn to fronts propagating into an unstable 
state for /j, > (for simplicity of notation, we use u = 0), and investigate the 
behavior of the integrand for £ — > oo. The large £ asymptotics of $o(0 follows 
directly from the o.d.e. fllTf ) by noting that /'($o(oo)) = /'(0) = fi, so that 



f»i j Ai(v) e~ A ~ € + A 2 {v) e" A+? , v > v* = , 



where 



X ± (v) = ^±y v i-4fi = l±ly/ v *-(v)2 , for//>0 (29) 

A*( V *) = A±(t;*) = y. (30) 

The behavior of <3>o for f = f * results from the fact that precisely at the 
so-called pulled velocity v*, the two roots A± coincide. 

While for an arbitrary velocity v > v* the term A\(v) in ( P8f) will be 
nonzero, so that the asymptotic behavior of $o is as e _A ~^, the pushed front 
solution — if it exists — is precisely the solution with a well-defined value 
v — v> at which Ai(i^) = 0. Note that for \i < 0, we have A_ < 0, so that 
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the relevant front solution in the range /i < has A\(v) = 0; thus the pushed 
front solution for \x > is precisely the analytic continuation of this front 
solution to the regime /x > 0. If such a solution with Ai(v^) = exists, it 
is the dynamically selected one from steep initial conditions |3l|. Moreover, 
these solutions decay for £ ^> 1 as e~ A +^, i.e., faster than e~ v As a result, 
integrands in (|7|) like (<%$ ) 2 or e v ^ (<%$ ) g(%^) oc e v ^ (<9 e $ ) $ for 
£ — > oo are integrable, as 

^t^^o^ €^ e ^t 5e _ 2 A +5 = e -VM) 2 -(«*) 2 C o. (31) 

Thus, for a pushed front both criteria for a solvability analysis of a pertur- 
bation theory are satisfied: the spectrum of the stability operator is gapped 
and the solvability integrals converge properly. 

In passing, we note that the adjoint mode x itself does not decay to zero 
for large £ in the supercritical range \x > 0, since 

€>i wt £ <9$ («t-V(«t) 2 -(«*) a ) C/2 5 ^oo 

v ~ e ? — — ~ e V / — ► oo . 32 

For our perturbation theory this is no problem as long as the inner product 
that defines the adjoint operator converges for £ — > ±oo. Eq. fl3~l|) shows that 
this is indeed the case. 

While the solvability integrals converge properly for pushed fronts, they do 
not for pulled fronts, as according to flSUp 

e ^(lf) 2 - 1 £ 2 e^e- 2 ^ = £ 2 ^oo. (33) 

As we already anticipated from the power law relaxation of pulled fronts, 
standard perturbation theory used to derive a MBA does not apply to pulled 
fronts. 

One could, of course, regularize the solvability integrals by first introducing 
a cutoff £ c , and taking the cutoff to infinity as the end of the calculation |63 



Whether such an approach yields sensible results, depends on the situation 
under consideration. If, e.g., this procedure is applied blindly to a solvability 
expression of the type (p7|) , one finds that the changes in the nonlinear terms 
of the equation give no contribution — in fact, since only the divergent terms 
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survive, this procedure amounts to calculating the changes in v* in perturba- 
tion theory for changes in the parameters in the linearized equation. Since v* 
can more easily be calculated explicitly from Eq. QTJ) such a calculation has 
no particular value. 

In fact, the divergence of the solvability integrals and the absence of a 
characteristic time scale for the internal front dynamics are deeply related. 
From (|25|) it is easily seen that the continuous spectrum defined by L <p a = 
—a 0O-, is bounded from below by a = [v 2 — v* 2 )/A. For a(k) = cu + k 2 , the 
eigenfunctions take the form of Fourier modes <fi a {k) e ±lk ^ in the leading edge 
region £ 3> 1. Hence for v = v*, the gap cro of the spectrum vanishes, and all 
the eigenfunctions of L are essentially plane waves in the semi-infinite leading 
edge. One finds furthermore that generic perturbations of pulled planar 
fronts $ are even outside the Hilbert space spanned by the eigenfunctions 
CT . In this case, the long time dynamics cannot easily be understood in terms 
of the eigenfunctions of L. One rather should directly study the linearized 
equation 



ed t 



e 2 V 2 + 



0=0 



+ 0(<j> 2 



(34) 



valid in the leading edge. In this formulation, the nonlinear region of the front 
interior plays the role of a boundary condition for the leading edge [EH], E7 



As a result one finds predictions like ([[]) - (f|). Note finally, that the leading 
edge extends on the same outer length scale on which also u varies. This 
demonstrates why it is not possible to eliminate the dynamics of a pulled 
front in a moving boundary approximation — independent of how thin the 
front is. 



5 Generalization of the solvability analysis and 
of its break- down 

In the previous sections, we have traced the main steps in the derivation of a 
MBA for two coupled equations that have been studied as phase field models 
for solidification. In this case, the inner equation for the order parameter 
reduces to the well-known nonlinear diffusion equation studied first by Fisher 
and Kolmogorov et al. 0, |65|, ^6|, and the nonhermitian linear operator L 



could be transformed to a hermitian operator L. This allowed us to obtain the 
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left zero mode x °f L explicitly. When one considers higher order dynamical 
equations or sets of coupled equations for the inner front region, it is usually 
not possible to find the adjoint mode explicitly Nevertheless, we show in this 
section that the same conclusions hold more generally. 

We consider where one has a vector <j>(x, t) of dynamical fields, 

that in the long time limit can approach a planar uniformly translating front 
profile $o(0 between the homogeneous stationary states ^ = $ (±oo). The 
front solution $o(0 with £ = x — vt obeys a set of o.d.e.'s, and because of 
translation invariance cf$o(£)/^£ is again a zero mode of the linear matrix 
operator L, obtained by linearizing the o.d.e.'s about the front solution $o(0 : 

If a front $o(0 is perturbed by external forces, other coherent structures 
or curvature effects, one generally encounters equations like 

L • 0i = gi (36) 

in a perturbation expansion about $o- in our example above, g\ decayed 
essentially like $ an d d§ /d£, as £ — > oo, and we only study such cases here. 
As is well known, such linear equations are solvable provided the right hand 
side is orthogonal to the kernel (null space) of L. The existence of a left zero 
mode x of L therefore generally leads to the solvability condition 

de*'& = , A^Q^, (37) 

(where the matrix Q contains some slowly varying fields), which relates pa- 
rameters of the expansion as in (|27|). So we now address the question of the 
existence of the left zero mode x °f L, which is defined through LJ • x — 0. In 
other words: x is the zero mode of the adjoint operator obtained by partial 
integration, 

/•OO POO 

d£b-{L-a)= d£ (L f • b) ■ a . (38) 



For this definition to hold, the integrals have to converge and the boundary 
terms that arise from performing the partial integrations all have to vanish. 
This imposes conditions on the allowed behavior of b, given the asymptotic 
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behavior of a: the product of these terms has to decay sufficiently rapidly for 

In general, there is no particular simplifying relation between L and L T ; 
e.g., a term f($ )d/d£ in L gives rise to a term —(df($o)/d£) — f($ )d/d£, 
in JJ. As a result, there is in general no simple relation between the left and 
right eigenmodes. However, since <&o(0 approaches the constant vectors (p^ 
for £ — > ±00, the operators L and asymptotically are linear operators with 
constant coefficients, so that 

4 f444-;*(f) , i = 4 / '-----^- 



^±00 « \d^ d^ 2 d^ 3 >rv\<*,j 



Moreover, in this limit, the operator L is exactly the same as the one that 
one obtains from linearizing the set of o.d.e.'s for <£> around the homogeneous 
stationary states (j)^. Therefore both $ an d the right zero mode d<&o/d£ of L 
are asymptotically for £ — > ±00 just sums of simple exponentials of the form 

d%/di X] ^ e ~ > (40) 

n=l 

where the eigenvalues Aj are determined by the characteristic polynomial of 
degree N 

det L( -A±, A± 2 , -A± 3 , ••• ; ± )=O. (41) 

The asymptotic behavior of an adjoint zero mode x follows immediately from 
the symmetry relation (|39|) . If we write the asymptotics for £ — > ±00 of x as 

X E^^- (42) 

71=1 

then the eigenvalues A^ are determined by the eigenvalue equation 

det V( - A± , A± 2 , -A± 3 , ■ • ■ ; $ ± ) = , (43) 
which in view of (|39| ) and the fact that det L = det immediately yields 

Aj = -Aj . (44) 
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Let us now investigate the asymptotic behavior of products b ■ a of left 
modes b and right modes a, which is required for the existence and definition 
of the adjoint operator and modes. Assume that the eigenvalues are ordered 
as Re A^ +1 > Re A^. A pushed or bistable front is a discrete solution with 
asymptotic behavior 

TV 

E 

=M- 
M 



A 



M+i e »*+ 



for £ 



oo 



n=M+l 
M 



(45) 



B M e- x M^ 



for £ 



-oo 



n=l 



that can be constructed from 



for a particular value of v 



;t. The 



right zero mode d<&o/d£ obviously has the same asymptotic behavior. In this 
expression, the eigenvalues Aj~, • ■ • , A^ f for £ — ► oo are all the eigenvalues with 
negative real part so that the exponentials converge, while on the right for 
£ — ► oo all A| f+1 , ■ ■ ■ ,X% have positive real parts. The existence of M modes 
on the left and N — M + 1 modes on the right is a reflection of the fact that 



the bistable or pushed front solution is an isolated (discrete) solution JH 



At this point, there is only one difference between bistable fronts and 
pushed fronts propagating into an unstable state: for the former, Re A < so 
that this mode is not present because it corresponds to a diverging behavior, 
while for a pushed front propagating into an unstable state, Re A| f > but 



by definition |3T 



A product of this right mode with a left mode converges to zero at ±oo, if 
the left zero mode behaves asymptotically like 



X 



M 



3 A n 5 



n=l 

N 



D M+1 e x M+it 



for £ 
for £ 



(46) 



-oo 



=M+1 



One easily verifies by counting the dimensions of stable and unstable manifolds 
in the two asymptotic regions, that also x belongs to a discrete spectrum, 
independent of the value of M, and that in general the divergent term ~ e A «^ 
is needed for this mode to exist. Indeed, the textbook argument 



L 



07 



Xl ■ 4>r. 



;l ! • xi) ■ a 



L f • Xi = 
Xi ■ (L • 



^ Xi 



Or. 



Xl 



(47) 
(48) 
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shows that the eigenvalues aj = a m equal each other, if the product of the 
eigenfunctions / xi ' 0m is finite and likewise that eigenfunctions with different 
eigenvalues are orthogonal. Application of simple "counting arguments" |56|, 
|3l| for the existence and multiplicity of solutions of o.d.eJs shows that (fj?D 
implies that associated with the discrete right zero mode of a pushed or bistable 
front solution, there is in general an isolated (discrete) left zero mode of 
with a nonzero divergent term (f4T|). 

This reasoning does not work for a pulled front, where the zero mode of L 
is part of a continuous spectrum with the same asymptotic decay properties at 
£ — > ±00. The same counting argument as above now yields, that in general 
no left zero mode of 1/ exists. 

This formal argument is supported by the observation, that a solvability 
integral for a pushed front diverges as the pushed velocity approaches the 
pulled velocity v* (|TJ): Generally, the velocity v will appear as a parameter in 
the characteristic polynomial (|4lD . If we consider the A n = A+ as functions 
of v, then according to the general scenario of front propagation into unstable 



states ]3T| the pulled velocity is associated with a minimum of the curve v(Xm) 
where Am is the root of ( f4T|) with the smallest positive real part. Hence for 
v ~ v* and uniformly translating fronts with Xm and Xm+i real, we have 



X M (v) = X* -y/v^V* -\ , (49) 

v 

2 



X M +i(v) = X* + —Vv-v* + ... , (50) 

V 

where 



v" 



d 2 v{X 



dX 2 M 



(51) 



is the curvature of v(Xm) in the minimum that determines v* and A* (see ||31|| , 
section V.C.2). It hence is a positive constant. 

In complete analogy with our earlier discussion in section 3, the general 
scenario for front propagation into unstable states is that while the asymptotic 
decay for £ ^> 1 is as e~ XkI ^ for an arbirary velocity v , a pushed front solution 
exists if for some velocity v' > v*, there is a front solution whose asymptotic 
large £ behavior is as e _Aju+1 ^ in agreement with (|45j). 

If there is no such pushed front solution, then starting from "steep" initial 
conditions the selected front velocity is v*; the asymptotic front profile with 
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this velocity is then 



$ (0 («e + /3)e- A «« , A A / = A A/+ i, (52) 

in analogy with (|28|). 

As we discussed above, for a pushed front, there is in general a discrete left 
zero mode fl46|) with asymptotic behavior \ ~ e AjW ^ for large £. In spite of this 
divergence, the product of left and right modes converges as 

X-Q-^| - 1 e^-^ C M • Q • A M+1 , 



~ e 



-(^'Ov^ 7 ? ^ o (v^v*) , (53) 



and solvability conditions generally can be derived. 

Just as we saw in the previous sections, the present analysis also shows 
that as v> approaches v* from above, the solvability integrals converge less 
and less fast until, at v* , we have according to (|52|) 



X ■ Q • - jr ~ C — > oo . (54) 

in complete analogy with our earlier result ([33|) for the example discussed in 
section 3. 



6 Conclusions and outlook 

In contrast to "bistable" or pushed fronts, the dynamics of pulled fronts is 
determined essentially in the leading edge. This was recently shown to imply 
a general power law relaxation of pulled fronts. In this paper, we have shown 
that this in turn entails that pulled fronts lack the separation of time scales 
necessary for the applicability of the usual MBA, and that solvability integrals 
diverge when a front is pulled. 

It is important to stress that one should not simply view this negative result 
as a formal problem — rather, one should take this conclusion as a signal that 
the pattern dynamics involving the motion of pulled fronts poses interesting 
new physical questions with possibly surprising non-standard answers. 
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As a first simple illustration of this, consider the uncoupled F-KPP equa- 
tion (0) in two dimensions with e = 1 and / = — 3 . If one starts with a 
radially symmetric steep initial condition, e.g., 4>(r,t = 0) = exp(— r 2 ), then 
this front will spread out in a circularly symmetric way. According to (|27|) the 
curvature correction will then give a contribution — 1/r = — l/(v*t) = — l/(2t) 
to the velocity at large times. However, in addition to that, there is a contri- 
bution —3/ (2b) of the same order of magnitude from the power law relaxation 
(0), as A* = 1 in this case. Thus, due to the combination of the power law re- 
laxation and the curvature correction, the front velocity v* will be approached 
asymptotically as v(t) = v* — 2/t p7| ! 

In the example above of a circularly symmetric pattern without any cou- 
pling to other fields, the relaxation and curvature effects can be simply added 
up, but for a less trivial patterns whose shape is changing in time, the proper 
description is far from obvious. Nontrivial patterns where pulled front propa- 
gation plays a dominant role occur e.g. in streamer discharges [[18]. Hence new 
analytical tools have to be developed for a moving boundary like description 
of these finger-like patterns. Work in progress suggests that the limit of 
zero electron diffusion creating shock-like electron fronts is a valuable approx- 
imation for negatively charged streamers. 

A recent illustration of the fact that the nonexistence of solvability integrals 
signals a transition to qualitatively different dynamical behavior is given by 



the behavior of fronts in the presence of multiplicative noise p9| , |30|| . Pushed 
fronts in the presence of multiplicative noise show regular diffusive behavior 
due to the noise being summed over the finite interior front region, and their 
diffusion coefficient can be expressed in terms of solvability type integrals [29 |. 
In contrast, fully relaxed pulled fronts in an infinite system do not diffuse at all, 
and if a front with pulled dynamics starts from a local (or "sufficiently steep" 
[ |31|j ) initial condition, it is subdiffusive PD[ : the root mean square displacement 
of pulled fronts increases with time as t 1 / 4 , not as t 1 / 2 . This prediction was 
first suggested by using a time- dependent cutoff £ c (t) ~ \ft in the solvability 
expression for the diffusion coefficient that is valid for pushed fronts. The 
motivation for this time-dependent cutoff comes from the relaxation analysis 
of pulled fronts given in Hence, this example illustrates both that using 



a cutoff in the solvability integrals sometimes can yield sensible results, and 
that the behavior of pulled fronts can be qualitatively different from those of 
pushed fronts. 

We finally note that these considerations also have implications for numer- 
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ical codes. In cases where a MBA applies in the limit in which the front width 
is taken to zero, numerical codes with adaptive gridsize refinement in the in- 
terior front region, where gradients are large, are quite efficient. For pulled 
fronts, however, solutions with a too coarse basic grid give inaccurate front 
velocities. For these, the refinement has to be done ahead of the front, in the 
leading edge 
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